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Abstract. We test a cosmological model which the only component is a pressureless 
5^ ' fluid with a constant bulk viscosity as an explanation for the present accelerated 

expansion of the universe. We classify all the possible scenarios for the universe 
I predicted by the model according to their past, present and future evolution and we 

' test its viability performing a Bayesian statistical analysis using the SCP "Union" data 

. set (307 SNe la), imposing the second law of thermodynamics on the dimensionless 

ly-j ' constant bulk viscous coefficient C and comparing the predicted age of the universe by 

CS| ■ the model with the constraints coming from the oldest globular clusters. 

I The best estimated values found for C, and the Hubble constant Hq are: C = 

^ ; 1.922 ±0.089 and Hq = 69.62 ±0.59 (km/s)Mpc"^ with a x^in = 314 (x^.^.f = 1.031). 

. The age of the universe is found to be 14.95 ± 0.42 Gyr. We see that the estimated 

(<~^ I value of Hq as well as oi Xd of very similar to those obtained from ACDM model 

^ ■ using the same SNe la data set. The estimated age of the universe is in agreement 

• i-H , with the constraints coming from the oldest globular clusters. Moreover, the estimated 

' value of ( is positive in agreement with the second law of thermodynamics (SLT). 

■ On the other hand, we perform different forms of marginalization over the 

parameter Hq in order to study the sensibility of the results to the way how Hq 
is marginalized. We found that it is almost negligible the dependence between the 
best estimated values of the free parameters of this model and the way how Hq is 
marginalized in the present work. 

Therefore, this simple model might be a viable candidate to explain the present 
acceleration in the expansion of the universe. 
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1. Introduction. 

In the last ten years the observations of type la supernovae (SNe la) have suggested 
a possible late time accelerated expansion of the universe (see for instance [l]-[7] and 
references therein). This discovery has additional support from the cosmic microwave 
background [8] and the large scale structure [9] observations. 

Several models have been proposed to explain this recent acceleration, one of them 
is the assumption of the existence of the so-called dark energy as responsible of such 
acceleration, i.e., a new unknown component that must be ~ 70% of the total content 
of matter and energy in the universe [l]-[7]. The leading dark energy candidates are 
a cosmological constant and a slowly varying rolling scalar field (quintessence models) 
[10]-[14]. However, a cosmological constant model faces several strong problems, one 
of them is the huge discrepancy between its predicted and observed value (of about 
120 orders of magnitude) [15]-[18], another one is the so-called the ^''cosmic coincidence 
problem" , i.e., why are we living in a moment where the matter density in the universe 
is of the same order than the dark energy density? [19] -[21]. 

On the other hand, it has been known since several years ago before of the discovery 
of the present acceleration in the expansion of the universe that a bulk viscous fluid can 
produce an accelerating cosmology (although it was originally proposed in the context 
of an inflationary period in the early universe) without the need of any cosmological 
constant or dark energy component [22]- [35] although some authors do not agree with 
this conclusion [36]. 

So, it is natural to think of the bulk viscous pressure as one of the possible 
mechanism that can accelerate the universe today (see for instance [3 7] -[46], where 
in some of these works the bulk viscous coefficient is assumed a priori without being 
derived from known physics or known particle properties). However, this idea faces the 
problem of that it is necessary to propose a viable mechanism for the origin of the bulk 
viscosity and in this sense some proposals have been already suggested [47]- [50]. 

Nowadays, with the observational data sets of SNe la [l]-[7] and some other 
cosmological observations like the shift parameter R of the cosmic microwave background 
radiation (CMB) given by the Wilkinson Microwave Anisotropy Probe (WMAP) 
observations [51], and the baryon acoustic oscillation (BAO) measurement from the 
Sloan Digital Sky Survey (SDSS) [52] it is possible to test cosmological models. So, 
in the present work we study and test bulk viscous matter-dominated cosmologies, 
i.e., scenarios which the only component of the universe is a pressurelcss fluid with 
bulk viscosity. The pressureless fluid characterizes to both baryon and dark matter 
components. The idea of this model is to explain the present acceleration of the universe 
using the bulk viscous pressure of the fluid without the need of any dark energy. 

As it was mentioned above, the explicit form of the bulk viscosity has to be assumed 
a priori or obtained from a known physical mechanism. In the present work we choose 
the first possibility. So, we assume a constant bulk viscous coefficient that it is the 
simplest parametrization for the bulk viscosity and we estimate its value using SNe la 
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observations. The theoretical solution of this model has been also analyzed in [22], 
[53]-[56] and other parametrizations for the bulk viscous coefficient C, have been also 
proposed and studied in [37]-[45] and [57]-[60]. There is also other general approach to 
the bulk viscous cosmologies called "fluids with inhomogeneous equation of state" (see 
[61, 65]). 

In section 2 we present the generalities of bulk viscous fluids in General Relativity 
(GR), in section 3 we apply this formalism to a bulk viscous matter-dominated universe 
model where we flnd the explicit expression for the Hubble parameter in function of 
the redshift. In section 4 we analyze the behavior of the scale factor for the possible 
scenarios that the model predicts for the universe according to the value of the bulk 
viscous coefficient. In section 5 we study the behavior of the deceleration parameter g, 
the curvature scalar i?, the total matter density pm and the age of the universe and in 
section 6 we briefly review the second law of thermodynamics (SLT). Section 7 presents 
the SNe la test to constrain the model and compute the best estimated values for the 
bulk viscous coefficient and the Hubble constant. After that, section 8 presents the 
best estimated values and the probability distribution functions for the bulk viscous 
coefficient by marginalizing over the Hubble constant assuming three different priors 
and finally in section 9 we present our conclusions. The appendices A, B and C detail 
the marginalization methods used in the present work. 

2. Theory of relativistic bulk viscous fluids. 

The origin of the bulk viscosity in a physical system is due to its deviations from the 
local thermodynamic equilibrium (for a review in the theory of relativistic dissipative 
fluids see [66]). In a cosmological fluid, the bulk viscosity arises when the fluid expands 
(or contracts) too fast so that the system does not have enough time to restore the 
local thermodynamic equilibrium and then it appears an effective pressure restoring 
the system to its local thermodynamic equilibrium. The bulk viscosity can be seen as 
a measurement of this effective pressure. When the fluid reaches again the thermal 
equilibrium then the bulk viscosity pressure ceases [49], [67]- [69]. 

In an accelerated expanding universe it is very possible that the expansion process 
is actually a collection of states out of thermal equilibrium in a small fraction of time. 
So, it is natural to assume the existence of a bulk viscous coefficient in a more realistic 
description of the accelerated universe today. 

We use the Weinberg formalism [70]-[72] for the imperfect fluids. So, in the present 
work we consider a bulk viscous fluid as source of matter in the Einstein flelds equations 
G^iv = STrGT^jy, where G is the Newton gravitational constant. 

On the other hand, the energy-momentum tensor of an imperfect fluid with a flrst- 
order deviation from the thermodynamic equilibrium has the form [49, 70, 71]: 

T^,^ = fm^u,, + {g^„ + u^u„)P* (1) 
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where 



P* = P- (yu, 



'V 



(2) 



In the equations (1) and (2) the four-velocity vector is of an observer who measures 
the effective pressure P*, P and p are the pressure and density of the fluid respectively. 
The term C is the bulk viscous coefficient that arises in the fluid which is out of the local 
thermodynamic equilibrium. 

It can be seen that the energy-momentum tensor (1) is similar to that of a perfect 
fluid but with an effective pressure P* composed by the usual pressure P of the fluid 
plus the pressure due to its bulk viscosity Pvisc = —CJ^^Uy. This viscous pressure, Pvisc, 
can be seen as a "measurement" of the pressure to restore the local thermodynamic 
equilibrium [49], [67]- [69]. The conservation equation for the viscous fluid is 



The effective pressure (2) was originally proposed by Eckart [74] in 1940 for a 
relativistic dissipative process in the context of thermodynamics systems out of local 
equilibrium, and subsequently Landau & Lifshitz presented an equivalent formulation 



However, the Eckart theory has problems in several aspects. One of them is that 
all the equilibrium states in this theory are unstable [76]. Another one is that signals 
can propagate through the fluids with superluminal velocities [77, 78] . 

In 1979, Israel- Stewart [79, 80] developed a more consistent and general theory that 
avoids these problems, and from which the Eckart theory is the flrst-order limit of the 
Israel-Stewart theory when the relaxation time goes to zero. Nevertheless, the Eckart 
theory is simpler to deal with than the Israel-Stewart theory. 

Despite of the inherent problems of the Eckart theory, and due to that it is simpler 
than the Israel-Stewart theory, it has been widely used recently by several authors to 
model bulk viscous dark fluids as responsible of the recent observed acceleration of the 
universe assuming that the approximation is valid for this purpose (see, for instance 
[37]-[41], [44, 45, 49, 56, 60]), i.e., in these papers assume a vanishing relaxation time, 
so that, in this limit the Eckart theory is a good approximation to the Israel-Stewart 
theory. In this context, it is convenient to mention that Hiscock et al [36] showed 
that flat Friedmann-Robcrtson- Walker cosmological models containing a bulk viscous 
Boltzmann gas expand more rapidly using the Eckart theory than the Israel-Stewart 
one, and suggesting that inflationary acceleration driven by bulk viscosity could be an 
effect of applying a pathological theory such as the Eckart theory. These results suggest 
that in this context the use of the Israel- Stewart causal theory would not produce a 
recent accelerating epoch as the Eckart theory would. However, posterior studies have 
suggested that this conclusion could not be true because consistent inflationary solutions 
have been found using the Israel- Stewart theory [34, 35]. 

It is important to point out that there exists other more general formulation for 
irreversible processes than the Israel-Stewart theory developed by D. Pavon et al. where 
the temperature is not necessary that of the thermodynamic equilibrium, see [81] for 



m"V,p+(p + P*)V"m, = 



(3) 



[75]. 
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details. In the present work we assume the Eckart theory and we constrain the model 
described above using SNe la observations. 

3. Cosmological model of bulk viscous matter-dominated universes. 

We study a cosmological model in a flat universe where the only component is 
a pressureless fluid with constant bulk viscosity as an explanation for the present 
accelerated expansion of the universe. The pressureless fluid characterizes both the 
baryon and dark matter components. 

Note that this model does not have the Cosmic Coincidence problem and that this 
fluid represents an unified description of the dark sector plus the baryon component in 
a similar way than the Chaplygin gas model (see for instance [44, 82] and references 
therein). In this approach the present acceleration of the universe is driven by the bulk 
viscous pressure of the fluid instead of a dark energy component. The present work 
differs from the approach of Cohstete et al [44] in that they propose two fluids, one of 
them is a bulk viscous fluid representing in an unifled way the dark sector with a bulk 
viscous coefficient C, proportional to a power of the energy density p (i.e., C — C i 
with C, 1/ constants), and the other one is a pressureless fiuid representing the baryon 
component. Nevertheless, most of their analysis and test using cosmological observations 
is done using the ansatz v = 0. 

Since we work with a pressureless fiuid (P = 0) then P* = Pyisc = —C^'^iJ'u, where 
( is the bulk viscous cocfiicient of the matter fiuid. 

We consider a spatially fiat geometry for the Priedmann-Robertson- Walker (FRW) 
cosmology as favored by WMAP [83] 

ds^ = -dt^ + (t) {dr^ + r^dn^) (4) 

where the function a{t) is the scale factor. On the other hand, we spread the conservation 
equation (3) in all its components 

Pm + (Pm - 3HC)3H = (5) 

where H = a/a is the Hubble parameter, pm is the total matter density, the dot means 
time derivative and W^u^, — 3H. 

The conservation equation (5) can be written in terms of the scale factor as 

= 3 {3HC - p^) (6) 

da 

This equation is valid for any parametrization of in particular we assume the ansatz 
C = constant that it is perhaps the simplest parametrization for the bulk viscous 
coefficient that can be proposed. Its value is going to be estimated from the SNe la 
observations. 

On the other hand, the flrst Priedmann equation for a flat universe is 

H = Pm (7) 
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So, substituting (7) into (6) we obtain 

dp, 



0, where 7 = 9 



c 



(8) 



da - ■■ ' ~ " V 3 ; 

Changing of variable from the scale factor to the redshift z using the relationship^ 
a — 1/(1 + z), we obtain the ordinary differential equation (ODE): 



-3Pm + 7P, 



1/2 







The exact solution of this ODE is: 



Pm(^) 



7 



+ 



{p]Ll - 1) (1 + ^) 



3/2 



(9) 



(10) 



where is the matter density evaluated today. Substituting this solution into (7) we 
obtain 

2 



H\z) 



Hi 



where Hq is the Hubble constant and where we have defined the dimensionless bulk 
viscous coefficient the matter density parameter and the critical density today 



Pcrit 



PmO 



Ho ~ ""^^ " 8nG 

In this model the bulk viscous matter is the only component of the universe implying 
that the first Fricdmann equation (7) evaluated today is Omo = 1- With this, the 
expression (11) finally becomes 



Pcrit — 



3H^o 



(12) 



= f [C+(3-C) (1 + ^) 



(13) 



4. Classification and evolution of bulk viscous matter-dominated models. 



In this section we classify all the possibilities for the universe predicted by this bulk 
viscous matter-dominated model using different values of the constant bulk viscous 
coefficient. Note that we have one different model for each value of so actually we 
have a collection of models depending of the value of (. 

We analyze the theoretical behavior of the scale factor in terms of the cosmic time. 
We begin expressing (13) in terms of the scale factor 



d Ho Ca'/^ + 3- C 



Integrating (14) yields 

Holt — to) — Ho 



r a"/^da' 
dt' ^3 / . 



(14) 



(15) 



where t labels the cosmic time and to the cosmic time today. 

I We assume that the value of the scale factor evaluated today is equal to one. 
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4.1. Case C = 0. 

When C = we recover the usual matter-dominated universe (with null bulk viscosity) 
with a scale factor coming from the integration of (15) like: 



a{t) = ( -i/o(i-io) + l 



2/3 



(16) 



The first and second derivatives of the scale factor with respect io x = HQ{t — to) are 



and 



da 
dx 

d'^a 



3a: + 2 



1/3 



9 1/3 



:i7) 



(18) 



The plot of functions (16)-(18) are shown in figures 1-3 respectively (the long 
dashed hues). 

This case predicts an eternal decelerated expanding universe (see figures 1-7). In 
this case the curvature scalar and the matter density are R = (3-?/g)a~^ and = 
{3HQ/87rG)a~^ respectively, that both diverge when the scale factor goes to zero (see 
sections 5.2 and 5.3 for details). When we have these two conditions (i?(a ^ 0) ^ 00 
and Pin (a ^ 0) ^ 00) in the past of the universe wc say that there was a Big-Bang. 
The elapsed time between the Big-Bang time till today is 

' (19) 



t. 



tn- 



where the subscript 'B' stands for "Big-Bang" . On the other hand, when a — > 00 then 
R and pm decrease to zero. 



4.2. Case C 7^0. 

For 7^ we do the change of variable y = a^/^ in expression (15) yielding: 



Ho{t - to) 



Cdy' 



This last expression can be rewritten as 

'c 



-.In 



Ca'/' + 3 - C 



3exp 



:^o(^ - ^o) 



ICV/^ + 3-cl 



(20) 



(21) 



When (a^^'^ -|- 3 — C > we can remove the absolute value bar of the right hand side 
term and then we obtain 

~-l 2/3 



ait) 



3exp[lCHo{t-to)) - 3 + C 



c 



Ct^O 



(22) 
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The first and second derivative of the scale factor with respect to x = i?o(^ — ^o) are 

-, 1/3 

(23) 



da 
dx 



C exp ( \Cx 



and 



d'^a 
dx'^ 



3exp(|Ca;) + C-3 



exp iUx] 2exp Ulx] + C - 3 



Sexpfka;) + C-3 



4/3 



(24) 



The behavior of the expressions (22)-(24) are shown in figures 1-3 respectively. 

The case Ca^/^ + 3 — C < does not correspond to a physical case of our interest 
because this case predicts an eternal contraction of the universe {H{t) < for any value 
of the cosmic time t) that is in contradiction with the observations (see the expression 
(14) to note that (a^/^ + 3 _ < q imphes H{t) < 0). 

In order to study carefully the expression (22) for different values of C, we take four 
cases: 

• < C < 3 

• C = 3 

• C>3 

• C<o 



4-2.1. Case < C < 3. From the expression (22) wc can sec that when t — > 00 then 
the scale factor tends to have the form like that of the de Sitter universe, i.e., 



a{t) oc e(^»^°(*-*°) 



(25) 



On the other hand, for any value of ( in this interval all the models predict an 
universe having a Big-Bang in the past in the cosmic time: 



tB — to 



(Ho 



In 1 



c 



(26) 



So, the universe begins with a Big-Bang followed by an eternal expansion (there 
is not any recoUapse epoch, see figures 1 and 2) and this expansion begins with a 
decelerated epoch followed by an eternal accelerated one (see figures 2 and 3). 

The transition between the decelerated-accelerated expansion epochs depends on 
the value of (. We compute the value of the scale factor where the transition happens. 
For that, we derive to d with respect to a using the expression (14): 



da 
da 




(27) 
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Figure 1. Plot of a{t,Q for dilferent values of Q in the interval < C ^ 3 (see 



expressions (16) and (22)). is the Hubble constant, t is the cosmic time and 
to is its present value. The long dashed line corresponds to C = (a flat matter- 
dominated universe with null bulk viscosity, see expression (16)). The short dashed 
line corresponds to C = 3 (the de Sitter universe). The blue line with the blue band 
correspond to a model with C, = 1.922 ± 0.089. This is the best estimated value of 
coming from the SCP "Union" SNe la data set analysis (see section 7). The band 
corresponds to the error at the 68.3% confidence level. 



Then, we make equal to zero to expression (27) in order to obtain the value of the scale 
factor "ot" where the transition happens, obtaining 



where the subscript "t" stands for "transition" . Writing the expression (28) in terms of 
the redshift z yields 



From the expression (28) we can see that for values of C, in the interval < C < 1 
the transition between the decelerated epoch to the accelerated one takes place in the 
future (at > 1). When — > 0, the value of at tends to infinity in the future. When 
= 1 then the transition takes place today (at = 1), when 1 < < 3 the transition 
takes place in the past of the universe (0 < at < 1) and when the value of Q is closer to 
3, the transition is closer to the Big-Bang time (see figures 2-4 and 7). 

The curvature scalar R is always positive and it has a very large value when the 
scale factor is very small (at the Big-Bang time) and it decreases forever tending to its 
minimum value R = |-f^oC when a ^ oo (see section 5.2 and figure 8) what is the 




(28) 




(29) 




Figure 2. Plot of the derivative of the scale factor with respect to x = Ho{t — to) for 
different values of C (see expression (23)). The long dashed line corresponds to C = 
(a flat matter-dominated universe with null bulk viscosity, see expression (17)) and the 
short dashed line corresponds to C = 3 (the de Sitter universe) . The blue line with the 
blue band correspond to a model with ( — 1.922 ± 0.089. This is the best estimated 
value of coming from the SCP "Union" SNe la data set analysis (see section 7). The 
band corresponds to the error at the 68.3% confidence level. 



de Sitter universe (consistent with the expression (25)). In the same way, the matter 
density pm is also very large when the scale factor is very small (at the Big-Bang time) 
and it decreases forever tending to its minimum value pm = {Hq/24:7iG)('^ when a — > oo 
(see section 5.3). 

The possible value for ( at least at the 99.9% confidence level estimated from 
the SCP "Union" SNe la data set is in the interval 1 < C < 3. So, assuming the 
best estimated value ( = 1.922 ± 0.089 (see table 1), the value of the scale factor at 
the transition time between deceleration-acceleration epochs is (see expression (28)): 
at = 0.42 ± 0.03 {zt = 1.33 ± 0.2). Note that the transition comes directly from the 
model, it is not necessary to assume any particular ansatz to estimate or induce it. 

4-2.2. Case C = 3. In the expression (22) we put = 3 to obtain 

a{t) = (30) 

This case corresponds to the de Sitter universe. Figures 1-5 show the behavior of 
expression (30) (the short-dashed lines). The scale factor becomes zero just when 
Ho{t — to) — oo and when Ho(t — to) oo the scale factor increases to infinity. 

In this case the model predicts an universe in an eternal accelerated expansion. 
The curvature scalar and the matter density are constants with values R = 12Hq and 
Pm = SHq/SttG respectively. See sections 5.2 and 5.3 for details. 
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X 



Figure 3. Plot of the second derivative of the scale factor with respect to a; = HQit—to) 
for different values of C, (see expression (24)). The long dashed line corresponds to C = 
(a flat matter-dominated universe with null bulk viscosity, see expression (18)) and the 
short dashed line corresponds to C = 3 (the de Sitter universe) . The blue line with the 
blue band correspond to a model with C, — 1.922 ± 0.089. This is the best estimated 
value of C coming from the SCP "Union" SNe la data set analysis (see section 7). The 
band corresponds to the error at the 68.3% confidence level. We see that for C < all 
the models predict a decelerated universe [d^a/dx^ < 0) forever. For C > 3 the models 
predict an accelerated universe {d^a/dx^ > 0) forever. And in the range < C < 3 the 
models predict a transition from deceleration to acceleration. 



Note that this model does not have a Big-Bang because the curvature scalar and 
the matter density are regular for any value of the cosmic time. 

4-2.3. Case ^ > 3. This case predicts an universe expanding forever, and this expansion 
is always accelerated (there is not any decelerating epoch or acceleration-deceleration 
transition). 

In this case when t ^ oo then the scale factor tends to the form of the de Sitter 
universe in the future (see expression (25)) and when t — >■ — oo then the universe tends 
to an Einstein static universe (defined by d, a = 0) where the value of the scale factor 
is not zero but tends to its minimum value: 



and during the evolution of the universe the scale factor is an increasing monotonic 
function with respect to the cosmic time. Thus, for this case there was not a Big-Bang 
and the age of the universe is not defined. Figure 5 shows some plots of the expression 
(22) for different values of > 3. 




) 



2/3 



(31) 
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Figure 4. Plot of i?(a, C) for different values of C (see expression (14)). a is the scale 
factor and Hq is the Hubble constant. The long dashed line corresponds to C = 
(a flat matter-dominated universe with null bulk viscosity) and the short dashed line 
corresponds to = 3 (the de Sitter universe). The blue line with the blue band 
correspond to a model with C, = 1.922 ± 0.089. This is the best estimated value of 
C, coming from the SCP "Union" SNe la data set analysis (see section 7). The band 
corresponds to the error at the 68.3% confidence level. 
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Figure 5. Plot of a{t, Q for different values of C, in the range C ^ 3 (see expression 
(22)). TJq is the Hubble constant, t is the cosmic time and is its present value. The 
short dashed line corresponds to C = 3 (the de Sitter universe). 

When the scale factor has the value a* the curvature scalar is zero and it is always 
increasing from zero until to reach its maximum value R = |-f^oC when a oo. See 
section 5.2 and figure 8 for details. 

The matter density pm has a similar behavior to the curvature scalar in the sense 
of that when the scale factor has the value a* the matter density is zero and it is always 
increasing from zero until to reach its maximum value pm = (-f^o/24vrG)C^ when a — > oo. 
See section 5.3 and figure 9 for details. 
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Figure 6. Plot of a(i, for different values of ^ in the range C < (see expressions 
(16) and (22)). Hq is the Hubble constant, t is the cosmic time and is its present 
value. The long dashed line corresponds to ^ = (a flat matter-dominated universe 
with null bulk viscosity, see expression (16)). 

4- 2.4- Case ( < 0. This case predicts an eternal decelerated expanding universe. The 
universe begins with a Big-Bang and expands forever until to reach its maximum value 
a* = (1 — when t — > oo, becoming an Einstein static universe in the future. The 

Big-Bang took place in the cosmic time defined by expression (26). Figures 2-4 and 
6 show the behavior of the scale factor for this case. 

In this case the curvature scalar has a transition from positive to negative values 
at = [(C — 3)/4C]^/^. When a — then i? — > oo (the Big-Bang) and when a = a* 
the curvature scalar is zero. The minimum value of the curvature scalar is negative and 
it is reached at the value of the scale factor a = [(2 — 6/C)/5]^/^. See section 5.2 and 
figure 8 for details. 

The matter density is always a decreasing function over its evolution starting from 
infinity (when a — > 0) to zero (when the scale factor is a*). See section 5.3 and figure 9 
for details. 



5. Cosmological parameters. 

5.1. Deceleration parameter q. 

We study the behavior of the deceleration parameter q function that is defined as 

, ^ aa d 1 , ^ 

The term d/a can be calculated from the second Friedmann equation, that for a matter- 
dominated universe with bulk viscosity reads: 

d AnG , , , , 

- = -^(Pm-9Ci^) (33) 
a 3 
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Prom the definition of C, (see expression in (12)) we have that 

On the other hand, from the first Friedmann equation (7) we have 

Pm = ^ H'{a) (35) 

So, substituting the expressions (34) and (35) in (33) we obtain: 
a 1 



^^CHo-H{a))H{a) (36) 



Thus, substituting (36) in (32) yields 

1 .Ho 



q{a) = - 1 - C 



2 V H{a) 
Using equation (14) for H{a) we arrive to 

Figure 7 shows the behavior of expression (37). Note that this expression docs not 
depend on Hq. We summarize the behavior of the deceleration parameter according to 
different values of ( as follows: 

• When ^ = we have q = \ that corresponds to a matter-dominated universe with 
null bulk viscosity. 

• When ( = 3 we have q = —1 that corresponds to de Sitter universe. 

• Case < ^ < 3: It is always a decreasing function from g(0) = | to q{oo) = — 1 
with a transition from positive to negative values in the value of the scale factor 
at = [(3 - C)/2Cf/^ (see expression (28)). 

• Case C > 3: When a — > a* then q — > —00 where a* — {1 — 3/C)^^^ is the minimum 
value of the scale factor (see expression (31)) and when a — > 00 then q — > —1. It 
is a negative and increasing function but it never becomes positive, its maximum 
value is -1. 

• Case C < 0: When a — then q — ^ and when a ^ a* then g — > 00 where 
a* — {1 — 3/C)^/^ is the maximum value of the scale factor. It is a positive and 
increasing monotonic function where its minimum value is |. 

From the expression (37) we obtain the deceleration parameter evaluated today as 

?(« = 1,C) = ^ (38) 

Note that when C = 1 the transition from the decelerated to accelerated epochs of 
the universe takes place today. For ( < 1 we have a decelerated universe in the present 
and for C > 1 we have an accelerated one today. 

Assuming the best estimated values for ( from table 1, the deceleration parameter 
today is (see expression (38)): gtoday = —0.46 ± 0.04 . 
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Figure 7. Plot of deceleration parameter q{a^ C) for different values of C (see expression 
(37)). The long dashed line corresponds to C = (a flat matter-dominated universe 
with null bulk viscosity). The short dashed line corresponds to C = 3 (the de 
Sitter universe). The blue line with the blue band correspond to a model with 
C — 1.922 ± 0.089. This is the best estimated value of ( coming from the SCP "Union" 
SNe la data set analysis (see section 7). The band corresponds to the error at the 
68.3% confidence level. 



5.2. The curvature scalar R. 



We calculate the curvature scalar R in order to study the singularities of the model. For 
a flat universe it is defined as 

a 



R = 6 



a 



Using the expression (36) in (39) we obtain 



R 



Hia) 



(Ho + 3H{a 

Substituting the expression (14) for H{a) into (40) yields 

■(3-0' , 5C(3-C) 



,3/2 



+ 4C 



(39) 



(40) 



(41) 



Figure 8 shows some plots of expression (41). From this expression we note that: 

• When C = then R = SH^/a^ that corresponds to a matter-dominated universe 
with null bulk viscosity. 

• When C = 3 then R = 12Hq that corresponds to de Sitter universe. 

• Case < C < 3: When a ^ then i? — > oo (the Big-Bang singularity). It is always 
a positive and decreasing function until to reach its minimum value R = |-f^oC 
when a oo. 

• Case ( > 3: It is always a positive and increasing function from zero at a = a* 
(where a* = (1 — 3/0^/^ is the minimum value of the scale factor) until to reach 



its maximum value R = ^Hq( when a ^ oo. 
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Figure 8. Plot of the curvature scalar i?(a, Q as a function of the scale factor for 
different values of ^ (see expression (41)). The long dashed line corresponds to ^ = 
(a flat matter-dominated universe with null bulk viscosity). The short dashed line 
corresponds to = 3 (the de Sitter universe). The blue line with the blue band 
correspond to a model with ^ — 1.922 ± 0.089. This is the best estimated value of 
C coming from the SCP "Union" SNe la data set analysis (see section 7). The band 
corresponds to the error at the 68.3% confidence level. 




Case C < 0- When a ^ then i? — > oo (the Big-Bang singularity). The curvature 
scalar is zero at the values of the scale factor ao and a*, where 

2/3 

ao = I — ^ I (42) 



4C 

and a* is the maximum value for the scale factor defined in expression (31). At 
there is a transition from positive to negative values of R. When < a < Oq the 
curvature scalar is a positive and decreasing function. For < a < a* the value of 
R is negative. The minimum value of the curvature scalar is i? = —^H'^C^'^ at the 
value of the scale factor a 

-| 2/3 

(43) 




5.3. The matter density p^. 

We calculate the matter density pm predicted for this model. We substitute the Hubble 
parameter H{a) (expression (14)) into the expression (35) for the matter density, yielding 

_n 2 

3-C 



Pm(a,C) 



Hi 
247rG 



c + 



,3/2 



(44) 
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Figure 9. Plot of the matter density pm(a,C) as a function of the scale factor for 
different values of C, (see expression (44)). The long dashed line corresponds to C = 
(a flat matter-dominated universe with null bulk viscosity). The short dashed line 
corresponds to C = 3 (the de Sitter universe). The blue line with the blue band 
correspond to a model with C, — 1.922 ± 0.089. This is the best estimated value of 
C coming from the SCP "Union" SNe la data set analysis (see section 7). The band 
corresponds to the error at the 68.3% confidence level. 

Figure 9 shows some plots of expression (44). From this expression we note that: 

• When C = then pm = (3//q / SttG) a^^ that corresponds to a matter-dominated 
universe with null bulk viscosity. 

• When C = 3 then = 3Hq/8ttG that corresponds to de Sitter universe. 

• Case < ( < 3: When a then oo (the Big-Bang singularity). It is 
a decreasing function until to reach its minimum value pm = {Hq/24:hG) ("^ when 
a — i> oo. 

• Case C > 3: When a = a* then pm(a*) = where a* = (1 - 3/C)^/^ is the 
minimum value of the scale factor (see expression (31) and section 4.2.3). It is 
an increasing monotonic function from zero until to reach its maximum value of 
Pm = {Hq/24:7^G) ("^ when a — » oo. In this case, the matter density increases as in 
the Hoyle's steady state cosmology [84]. 

• Case C < 0: When a — > then pm ^ oo (the Big-Bang singularity). It is a 
decreasing function until to reach its minimum value in pm(fl*) = 0. 
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5.4.. The age of the universe. 

We define the age of the universe as the elapsed time between the Big-Bang time 
until today to. So, from the expression (19) for the case C = we have 

Age = |tB - tol = for C = (45) 

And from the expression (26) for the case C < 3 

Age= I^B-tol = -~^ln I 1- I I , for C<3 (46) 



Note that for the case C, > ^ the age of the universe is not defined. On the other 
hand, evaluating the expression (46) in the best estimated values of C and Hq from the 
SCP "Union" SNe la data set (see table 1) yields an age of 14.957 ± 0.422 Gyr. The 
errors are at 68.3% confidence level. We found that this age is in perfect agreement with 
the constraints on the age of the universe coming from the oldest globular clusters [85] . 
Figure 10 shows the age of the universe for different values of C, and the best estimate, 
the vertical hues correspond to Hq — [55, 75] (km/s)Mpc~^ §, it is the permitted region 
according to values of Hq consistent with the distance moduli used to derive ages for 
Galactic globular clusters from the Hipparcos parallaxes. As a reference, the predicted 
age of the universe for ACDM (Lambda Cold Dark Matter) cosmological model with 
a fiat universe and using the best estimated values of flmo and Hq shown in table 1 is 
13.750 ± 0.29 Gyr. 

6. Thermodynamics and the local entropy. 

The law of generation of the local entropy in the FRW space-time is found to be 
[70, 71, 72] 

TV^^s"" ^CV^u'' ^3HC (47) 

where T is the temperature and V,^s^ is the rate at which entropy is being generated in 
an unit volume. Then, the second law of the thermodynamics can be written as 

TV^s" > (48) 

which, from the expression (47), it implies that 3H( > 0. 

Since the Hubble parameter H is positive in an expanding universe then C has to 
be positive in order to preserve the validity of the second law of the thermodynamics. 
Thus, equation (48) can be written for this model as 

C > (49) 

§ Hereinafter the units of Hq are expressed in (km/s)Mpc~^ 
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Figure 10. Plot of the age of the universe in units of Gigayears (Gyr) with respect 
to iJo in units of (km/s)Mpc^^, for different values of C, (see expression (45), (46)). 
The blue point is located at 14. 957± 0.422 Gyr that corresponds to the best estimated 
value for the age of the universe coming from the SCP "Union" SNe la data set (where 
Hq = 69.62 ± 0.59(km/s)Mpc~^, see section 7). The blue line with the blue band 
correspond to a model evaluated at the best estimated value for C, = 1.922 ± 0.089. 
The vertical lines correspond to the interval Hq = [55, 75] (km/s)Mpc~"^, it is the 
permitted region according to values of Hq consistent with the distance moduli used 
to derive ages for Galactic globular clusters from the Hipparcos parallaxes (see ref. 
[85]). The horizontal lines corresponds to the constraint for the age of the universe 
from the oldest globular clusters (Age= 12.9±2.9 Gyr [85]). So, the shaded area is the 
consistent region for the age of the universe. The dashed line corresponds to C = 0, a 
flat matter-dominated universe with null bulk viscosity. 



7. Type la Supernovae test. 

We analyze and constrain the viability of the model using the type la supernovae (SNe 
la) observations and the entropy tests (see expression (49) and section 6). For that, 
we calculate the best estimated values for the parameters ( and Hq and the goodness- 
of-fit of the model to the data by x^-minimization and then compute the confidence 
intervals for {(, Ho) and the probability density functions for ^ (marginalizing over Ho) 
to constrain their possible values. 

We perform the statistical analysis using the "Union" SNe la data set of "The 
Supernova Cosmology Project" (SCP) [73] composed by 307 type la supernovae brought 
together from 13 independent data sets. 

We use the definition of luminosity distance [3, 5], [86]-[88] in a fiat cosmology, 

it is 




(50) 



Constraining a matter- dominated cosmological model with constant bulk viscous 



20 



10 



+ 25 (51) 



where H[z, Hq) is the Hubble parameter (expression (13)) and 'c' is the speed of light. 
The theoretical distance moduli for the k-ih. supernova with redshift is defined as 

Mpc 

where m and M are the apparent and absolute magnitudes of the SNe la respectively 
and the superscript 't' stands for theoretical. We construct the statistical function as 

r ~ 1 ^ 

^ \li\zkX,Ho) - Hk 
x'iC, Ho) = -2 (52) 

where /ik is the observational distance moduli for the k-th supernova, is the variance 
of the measurement and n is the amount of supernovae in the data set. 

Once constructed the function (52), we numerically minimize it to compute the 
"best estimates'^ for the free parameters of the model (i.e., (,Hq). The function 
measures the goodness-of-fit of the model to data. The probability density function 
(pdf) is defined as 

pdf(C, Ho) = cte ■ e-^'/2 ^53) 

where 'cte' is a normalization constant. 

So, we minimize the function to obtain the results shown in table 1. The 
confidence intervals for ( and Ho are shown in figure 11. Note that at least at the 
99.73% confidence level the possible values for ( are inside of the interval 1 < ^ < 3. 
The models in this interval indicate the existence of a Big-Bang in the past of the 
universe, they have an accelerated expansion today and predict an age of the universe 
~ 13 — 18 Gyr at the 99.73 % confidence level (see figure 11 and subsections 4.2.1, 5.4). 

In order to compare the results shown in table 1 (in particular the measurement 
of the goodness-of-fit to data Xd o f ^^^^ host estimated value for Ho) with other 
models we compute also the best estimated values for the ACDM model using the same 
SNe la data set. For the ACDM model, wc assume a fiat cosmology with Ho and 
r2i„o as free parameters (implying that the cosmological constant density parameter is 
f^Ao = 1 — ^mo)- The results are shown in Table 1. 

We find that the best estimated values for Ho and the goodness-of-fit to data x^.o-f. 
from ACDM model are very similar to those obtained from the bulk viscous model using 
the same SCP "Union" SNe la data set and in agreement with that reported by 5 year 
WMAP {Ho = 70.5 ± 1.3 (km/s)Mpc-\ see [83]). 



8. Marginalization over Ho 

Now, we construct a pdf for the bulk viscous coefficient C marginalizing over the Hubble 
constant Ho in order to have to C as the only free parameter of the model. We use three 
different priors to marginalize Ho'. 

(i) Constant prior over Ho (see Appendix A). 



Constraining a matter- dominated cosmological model with constant hulk viscous 21 

Confidence intervals 




Figure 11. Confidence intervals for C and Hq of a bulk viscous matter-dominated 
universe model with constant bulk viscosity. The best estimated values and confidence 
intervals were computed using the SCP "Union" SNe la data set where it is found 
C = 1.922 ±0.089 and Hq = 69.62 ±0.59 (km/s)Mpc~^ as best estimates (see table 1). 
The shown confidence intervals correspond to 68.3%, 95.4% and 99.73% of probability 
for the estimated value. Hq is in units of (km/s)Mpc~"'^ and C is dimensionless. Note 
that at least at the 99.73% confidence level the possible values for C are inside of 
the interval 1 < C < 3 indicating an accelerating universe today with high confidence 
level (see subsection 4.2.1). The blue dashed lines show different ages of the universe 
according to the values of C, and Hq (see expression (46)). The red vertical dashed 
lines indicate values of the deceleration parameter evaluated today (see expression 
(38)) being the line go = the limit between a decelerated-accelerated universe today. 
The estimated age of the universe is 14.957 ± 0.422 Gyr, where the error is at 68.3% 
confidence level. We found that this age is in perfect agreement with the constraints 
on the age of the universe coming from the oldest globular clusters [85] (see figure 10). 
The models that are in the grey region {C, > 3) do not have a Big-Bang in the past of 
the universe as the label indicates. 



(ii) Gaussian prior over Hq centered at Hq = 70.5 ± 1.3 (km/s)Mpc ^ and Hq = 
72 ± 8 (km/s)Mpc~^ as reported by 5 year WMAP data [83] and HST Cepheid 
variable star observations [89] respectively (see Appendix B). 

(iii) Dirac delta prior over Hq centered at Hq = 70.5 (km/s)Mpc~^ and Hq = 
72 (km/s)Mpc~^ reported as mentioned above (see Appendix C). 



8.1. Marginalization assuming a constant prior over Hq. 

We marginalize over Hq by assuming a constant prior as described in Appendix A. 
We use the Xcp function (A. 18) instead of function (52) to perform the supernova 
analysis. The results are: 

. C = 1.922 ±0.089, with a xLn = 331.30 (xL.f. = 1.08). 

The error is at 68.3% confidence level. We plot the pdf of ( to show the value that 
maximizes such a function, i.e., the best estimated value of This result is shown in 
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Viscous model C, — constant 



Model 


Ho 


c 




A mm 


Xd.o.f. 


Constant bulk viscosity 


69.62 ±0.59 


1.922 ± 0.089 


1 


314.57 


1.031 


ACDM 


70.01 ± 0.59 




0.278 ±0.027 


311.84 


1.022 



Table 1. Summary of the best estimates of C, and Hq for the bulk viscous matter- 
dominated universe model with a constant bulk viscosity coefficient. Also, it is shown 
the best estimates for the ACDM model in order to compare the values of Hq and Xmin 
of both models. It is found that the best estimates for Hq and their corresponding x^j^^ 
computed for both models and using the same SNe la data set are almost the same 
value and in agreement with the reported value by 5 year WMAP {Hq = 70.5 ± 1.3 
(km/s)Mpc^^, see [83]). For the ACDM model, we assumed a flat cosmology with 
Hq and fimo as free parameters, where Om is the matter density parameter. Note 
that in the constant bulk viscous model the value of ri,„ is not estimated but it is 
assumed to have the value of one as part of the model. The best estimates were 
computed by a Bayesian statistical analysis using the SCP "Union" 2008 compilation 
data set composed by 307 SNe la from 13 independent data sets [73]. iJo is in units of 
(km/s)Mpc~^ and C,, f^m are dimensionless. The subscript "d.o.f." stands for degrees 
of freedom and the errors are at 68.3% confidence level. Figure 11 shows the confidence 
intervals. 



Viscous model C = constant 

Best estimates for C, 



Marg. 


c 


Amin 


Xd.o.f. 


(i) 


1.922 ± 0.089 


331.30 


1.08 


(ii) 


1.941 ± 0.084 


315.26 


1.03 


(iii) 


1.924 ±0.089 


314.91 


1.02 


(iv) 


2.026 ±0.051 


316.75 


1.03 


(v) 


2.195 ± 0.049 


330.41 


1.07 



Table 2. Summary of the best estimates of C for the bulk viscous matter-dominated 
universe model with a constant bulk viscosity coefficient. The best estimates were 
computed by a Bayesian statistical analysis using the SCP "Union" 2008 compilation 
data set composed by 307 SNe la from 13 independent data sets [73]. The errors are 
at 68.3% confidence level and ( is dimensionless. Figure 12 shows the probability 
density functions for (. We have marginalized over the Hubble constant Hq using 
three different priors described in section 8. In the first column of the table the label 
"Marg." indicates the marginalization method over Hq used to obtain the results. 
The numbers correspond to: 

{i) Constant prior over Hq (see Appendix A). 

{ii)Gaussian prior over Hq centered at Hq = 70.5 ± 1.3 (km/s)Mpc~^, as reported 

by 5 year WMAP data [83] (see Appendix B). 
(iii) Gaussian prior over Hq centered at -ffo = 72 ± 8 (km/s)Mpc~^, as reported by 

HST Cepheid variable star observations [89] (see Appendix B). 
{iy)Dirac delta prior over Hq located at Hq = 70.5 (km/s)Mpc~^, as reported by 5 

year WMAP data (see Appendix C). 
{v)Dirac delta prior over Hq located at Hq = 72 (km/s)Mpc~^, as reported by HST 
Cepheid variable star observations (see Appendix C). 
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table 2 and figure 12. 

8.2. Marginalization assuming a Gaussian prior over Hq. 

We marginalize over Hq by assuming a Gaussian prior as described in Appendix 
B. We use two different central values for the Gaussian prior for Hq, they are 
Hq = 70.5 ± 1.3 (km/s)Mpc~^ and Hq = 72 ± 8 (km/s)Mpc~^ as suggested by 5 year 
WMAP data [83] and HST Cepheid variable star observations [89] respectively. See 
section Appendix B for details. The results are shown in table 2 and figure 12. 

8. 3. Marginalization assuming a Dirac delta prior over Hq . 

Finally, we marginalize over Hq by assuming a Dirac delta prior as described in 
Appendix C. In practice, it means to assume a specific value for Hq. In the same 
way as above we use two different values for Hq, they are Hq = 70.5 (km/s)Mpc~"'^ and 
Hq = 72 (km/s)Mpc"^ The results are shown in table 2 and figure 12. 

9. Conclusions. 

We performed a detailed study of a bulk viscous matter-dominated universe model 
with constant bulk viscous coefficient (. The only component of this model is a 
pressureless fiuid with an assumed constant bulk viscosity driving the present accelerated 
expansion of the universe. We analyzed the different possible scenarios for the behavior 
and evolution of the universe predicted by the model according to the value of the 
dimensionless bulk viscous coefficient C. In general, the model predicts an expanding 
universe for any value of C and we find that for values of C in the range < C < 3 
the model predicts an universe with a Big-Bang in the past and then starting with a 
decelerated expansion epoch followed by a transition to an accelerated epoch in late 
times. 

The case C = 3 corresponds to the de Sitter universe. When C > 3 the universe is 
always in an accelerated expansion and the value of the scale factor tends to infinity in 
the future (in a similar way than the de Sitter universe) and where the curvature scalar 
goes to a positive constant value. In this case there is not a Big-Bang nor a defined age 
nor origin of the universe because the scale factor has a minimum value greater than 
zero in the past of the universe (when t — > — oo), where also, the curvature scalar and 
the matter density are zero, i.e., we have the Einstein static universe in the past. 

When C < the universe is always in a decelerated expansion and the scale factor 
tends to a finite value in the future (it tends to the Einstein static universe in the future) 
where also, the curvature scalar and the matter density go to zero. For C, in this range 
there is always a Big-Bang in the past. 

For any value of C there does not exist any recoUapse epoch or transition between 
expansion to contraction epochs. When C = 1 the transition between a decelerated to 
accelerated epoch takes place today (i.e., qo — Oioic( — l). For < C < 1 the transition 
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Figure 12. Probability distribution functions for the coefficient Q in 8. bulk viscous 
matter-dominated universe model with a constant bulk viscosity. We have marginalized 
over the Hubble constant Hq using three different priors described in section 8. The 
best estimates and pdf's were computed by a Bayesian statistical analysis using the 
SCP "Union" 2008 compilation data set composed by 307 SNe la from 13 independent 
data sets [73]. Table 2 summarizes the best estimated values of C, using the different 
marginalizations. The first figure (on top) corresponds to the marginalization using a 
constant prior over Ho (see Appendix A). The second and third figures (from top to 
bottom) correspond to the marginalization using a Gaussian prior over Hq centered 
at Ho = 70.5 ±1.3 (km/s)Mpc~^ and iJo = 72 ± 8 (km/s)Mpc~^ as reported by 5 year 
WMAP data [83] and HST Cepheid variable star observations [89] respectively. The 
fourth and fifth figures correspond to the marginalization using a Dirac delta prior 
over Ho located at Ho = 70.5 (km/s)Mpc~^ and Hq — 72 (km/s)Mpc~^, as mentioned 
above. 
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between the deceleration-acceleration epoch takes place in the future and for 1 < < 3 
the transition took place in the past. 

On the other hand, we estimate the values of C and using the SCP Union SNe 
la data set where it is found that the estimated values of C, are positives, in agreement 
with the second law of thermodynamics. We note that the best estimated value for Hq 
is in agreement with that reported by 5 year WMAP and estimated from ACDM model 
(using the same SNe la data set) and the minimum value of the function {x\.o.i) 
obtained is of the same order than that obtained for ACDM model indicating that this 
is a competitive model to fit the supernova observations. 

We estimate and compare the age of the universe predicted by this model with the 
constraints on the age coming from the oldest globular cluster observations. We find 
that the best estimated values for the age of the universe is 14.95 ± 0.42 Gyr that is in 
agreement with the constraint in the age coming from the oldest globular clusters. 

The present model has been constrained using SNe la data. It means that we 
have constrained the model using information (observations) of the past of the universe 
until a redshift z < 1.5. So, SNe la observations are not enough to trace the whole 
evolution of the universe from the Big-Bang until the present time. Therefore, it would 
be important to constrain the model using also cosmological observations that can 
provide information of the early universe like the observations of the CMB anisotropics 
from the WMAP experiment and the large scale structure (LSS) from SDSS experiment. 
A preliminary work in that sense has been already done in [56] showing (but not in a 
conclusive way) that when the model is constrained using the shift parameter R of 
CMB and the parameter A measuring the BAG peak from SDSS, in addition to the 
Gold 2006 SNe la data set, the estimated value of C, is negative at 99.7% confidence 
level, violating the second law of thermodynamics and the estimated value of Hq is 
lower (ifo ~ 53(km/s)Mpc"^) than that reported by 5 year WMAP and HST Cepheid 
variable star observations and giving a bad goodness-of-fit to data. Considering these 
facts we may conjecture that this model does not work well for early times of the universe 
and that in order to be a viable model the bulk viscosity should be triggered just until 
late times. 

Another problem of the model is to explain the origin of the bulk viscosity from 
known or new physics. Some proposals in that sense have been explored, for instance 
in [49] and [50] , where it is proposed a mechanism to generate the bulk viscosity by the 
decay of dark matter particles into relativistic products. 

Finally, we explore also the sensitivity of the results to different ways of 
marginalization over when we estimate just the coefficient C,. We find that it is 
almost negligible the difference in the estimations of C, when Hq is marginalized assuming 
constant and Gaussian priors centered in the values reported by 5 year WMAP and HST 
Cepheid variable star observations. However, when Hq is marginalized by assuming a 
specific value (a Dirac delta prior located at the values reported as mentioned above) 
the estimated values of C have a small increase and the values of x\o.i. increase as well 
(i.e., we have a worse fit to data). 
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Appendix A. Marginalization assuming a constant prior over Hq. 

To construct a pdf that depends only on the parameter C, we use the process of 
marginalization over Hq in order to eliminate the dependence of the pdf with respect 
to the parameter Hq. 

The following procedure can be applied for any other model. It is, in a model with 
several parameters where it is necessary to reduce the number of free parameters or, 
for some reason, to eliminate the dcpendancc of the model from some particular free 
parameters (for instance: Hq) in order to compute the best estimated values of the other 
free parameters that we keep (in our case: Q the solution is to marginalize over the 
parameters that we want to eliminate. 

In the present work we assume three different prior distribution functions for Hq: 
constant, Gaussian and Dirac delta. In this appendix we describe the marginalization 
using a constant prior for Hq. 

We start by multiplying and dividing by Hq the equation (50) 

ddzX.HQ)^^—^ ; (A.l) 

-Ho Jo H{z',C,Ho) 

Next, we define a new dimensionless luminosity distance as Dl{z, Q = HQ-dL{z, C, -f^o) / c. 
Then 

D,{zX)=c{l + z) (A.2) 
Jo E{z',Q 

where Dl{z, does not depend on Hq anymore and E{z, = H{z, C, Hq)/ Hq. 
The theoretical distance moduli (51) becomes 

t,\z, C, Hq) = 5 log,o (^ ^hI'mpc ) + ^ ^ (^^) ^ ^^-^^ 

where we have defined a dimensionless "Hubble parameter" Hq = Hq ■ Mpc/c. It is 
useful to define a new theoretical distance moduli that does not depend on Hq anymore 
as 

~^\z,C)=5\og,Q[DL{z,C)] + 25 (A.4) 
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Then, the expression (A. 3) for the distance moduh becomes 

fx\z,C,Ho) = fl\z,C) - 51ogio[^o] (A.5) 
Now we construct the function (52) with these new definitions as 



2 



(A.6) 



i=l 



where lif^ is the observed distance moduh and CTj its variance. 
We rewrite the expression (A.6) as 



obs \ ^ / X " / r, t , , obs 

f-j H'i 



X^C ^o) = E ^^^-^ ) - 2 (5 log.o ^o) E 

i=l ^ * ^ i=l 



a? 



2 " / 1 \ 



+ ( 5 logio Ho 

If we define || 

/ ,~,t ..ODSX" '- ,~,t ,,ODS "- 1 



/ ,~,t ,,obs\ 2 n ~t ,,obs 



i=l ^ i=l ^« 1=1 

then we can express (A. 7) as 

X^ {C,Ho)^A- 2Bx + Cx'^ (A.9) 

where 

x = 51ogio(^o). (A.IO) 

Note that all the dependence of the function with respect to Hq is now in the x 
variable. Remind that the pdf is defined as 

pdf(C,^o) = cie-e-^'/' (A.ll) 

where "cie" is a normalization constant. So, we marginalize the pdf (A.ll) over Hq 
computing the following integration 

/CO 
pdi{CHo)-pdi{Ho)dHo (A.12) 
-oo 

where pdf(ifo) is the prior probability density function for Hq. In general, the 
integration (A.12) has to be done over the range of all the possible values of the 
parameter to be marginalized. In this case the range of possible values for Hq, or 
equivalently Hq, is (—00,00). 

If we take the case where the prior for Hq is a constant then: 

/oo ^ 
pdf(C,^o)c?^o (A.13) 
-00 

11 Note that these expressions do not depend on Hq anymore. 
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A constant prior means that we do not prefer any particular value for Hq, i.e., any 
value for Hq has the same probability of being. To solve the integral (A. 13) wc do a 
change of variable from Hq to x (see expression (A. 10)). At the same time, we substitute 
the expression (A. 9) into (A. 11) and then (A. 11) into (A. 13), yielding 



,'lnlO' 

pdf(C) = cte ( J exp 



1 



-A 



exp 



X — 



B 
C 



dx 



(A. 14) 



where B = B + (lnl0)/5. 

We can see that the integral (A. 14) has the form of a Gaussian distribution with 
value: 

{x — xY' 



— f 



exp 



oo 

2 



2(7^ 



dx 



(A.15) 



where x is the mean and a is the variance of the distribution. 
Therefore the expression (A. 14) becomes 

/lnlO\ [2^ 
pdf (C) = cte l—^ j V 




(A.16) 



so that the pdf ((^) does not depend on Hq anymore. Note that it was not necessary any 
numerical integration of the expression (A. 13). We can express this pdf (C) in terms of 
a new xlp function like 



pdf (C) = a ■ e-^'p/2 
where a = cte ■ V2nla 10/ (by/C) and 

, ^ ^ [5(C)+ln(10)/5'' 



(A.17) 



C 



(A. 18) 



This new xlp function does not depend on Hq anymore. The label "cp" stands for 
constant prior for Hq. 



Appendix B. Marginalization assuming a Gaussian prior over Hq. 



We perform the marginalization assuming a Gaussian probability distribution function 
for Hq centered at Hq and with standard deviation a*. So, the pdf (Hq) prior with the 
form of a Gaussian distribution for Hq is 

2" 



pdf(i/o)=exp 
With this, the expression (A. 12) becomes 

POO 

pdf(C) — cte ■ / e"-^ exp 



fB.ll 



I f Hn — Hn 



a' 



dHQ 



(B.2) 
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where the function is given by (A. 9) and "cte" is a normahzation constant. 

In the present work we use two different central values for the Gaussian prior 
coming from two different observations. One of them, is that coming from the 5 year 
WMAP observations where = 70.5 (km/s)Mpc~"'^ with an standard deviation of 
a* = 1.3 (km/s)Mpc^^. The other central value is Hq = 72 (km/s)Mpc~^ with an 
standard deviation of a* — 8(km/s)Mpc~^ as reported by the HST Cepheid variable 
star observations. 

In practice we perform the numerical integration of the expression (B.2) in the 
interval Hq — [55, 85] (km/s)Mpc~^ considering that this is a suitable and representative 
interval that include almost the 100% of the probability density for Hq. 

Appendix C. Marginalization assuming a Dirac delta prior over Hq. 

It means to assume a specific value of Hq. It is hke to think that the value of the 
Hubble constant is Hq that has been measured with an infinite accuracy (with standard 
deviation equal to zero), so that its pdf(ifo) has the form of a Dirac delta function. 

Clearly, this assumption is just an idealization to simplify the work and it does not 
correspond to the reality. However, for some cases it is a good approximation. This 
prior has the advantage of that once assumed it is very simple to perform the integration 
of the expression (A. 12) using the Dirac delta properties. Thus, the prior with the form 
of a Dirac delta for Hq is 



where "cte" is a normalization constant. 
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